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Abstract 

In the strong Yukawa coupling limit the renormalization of trilinear couplings Ai{t) 
and combinations of the scalar particle masses 971? (t) are studied in the framework 
of the Next-to-Minimal SUSY Model. The dependence of given parameters on 
their initial values at the Grand Unification scale disappears when the solutions 
of the renormalization group equations approach the infrared quasi fixed points. 
Besides that at dcuT ^ n(0) ^ 1 in the vicinity of quasi fixed points Ai{t) and 
971? (t) concentrated near lines or surfaces in the soft SUSY breaking parameter 
space. It occurs because the dependence of trilinear couplings and combinations of 
scalar particle mass on Aj(0) and 971? (0) disappears rather weakly. We propose a 
method which enables to obtain equations for such lines and surfaces. Their origin 
is also discussed by means of exact and approximate solutions of the NMSSM 
renormalization group equations. 
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1 Introduction 


Parameters of a soft breaking of supersymmetry plays a key role in analysing the particle 
spectra in supersymmetric (SUSY) models. The Bose-Fermi degeneracy of the spec¬ 
trum is one of the most serious flaws in SUSY models. This means that the masses of 
observable particles and their superpartners coincide in the limit of exact supersymme¬ 
try, but this is in glaring contradiction with experimental data. Thus, supersymmetry 
must be broken, but its breaking must not lead to the hierarchy problem [|^. Such a 
breaking of supersymmetry is referred to as a soft breaking. Supersymmetry breaking 
associated with supergravity (SUGRA) effects is one of the most promising mechanisms 
for constructing a realistic model. Although the Lagrangian of Y = 1 SUGRA models 
1^ is unrenormalizable, it can be shown that, in the low-energy region E -C Mpi, where 
Mpi = 2.4 ■ 10^® GeV is the Planck mass, all unrenormalizable terms in the Lagrangian 
for observable helds are suppressed in a power-law way with respect to Mpi, vanishing 
in the limit Mpi —> oo. In the case being considered, the Lagrangian for observable fields 
can be represented as the sum 


^ — -^SUSY + Aoft, 


( 1 ) 


where £susy is the Lagrangian corresponding to unbroken supersymmetry, while £soft 
takes into account terms that generate a soft breaking of supersymmetry. 

To a considerable extent, the form of the hrst term in (|^ is determined by the super¬ 
potential of the SUSY model under study, this superpotential in turn being a function 
of chiral superhelds. For a SUSY theory to be renormalizable, the superpotential must 
include only terms that are quadratic and cubic in chiral superhelds Sa- 

W{Sa) = -fJ-a/dSaS/B + -ha/B'ySaSpS^. ( 2 ) 

When local supersymmetry is broken in SUGRA models, parameters of a soft breaking of 
supersymmetry are generated in the sector of observable helds |§],0. Since the explicit 
form of terms that violate supersymmetry, but which do not lead to the emergence of 
second-order divergences, is known |P|, the Lagrangian for the case where supersymmetry 
is softly broken can be represented as 




mlylva 


Q-^afB'ihoifB'y'ya'yfBy'f T 2h.c. 


(3) 


where y^ are the scalar components of the chiral superhelds Sa and are gaugino helds. 
Usually, the parameters Mq, m^, and Bap of a soft breaking of supersymmetry 

are dehned at the scale of Mx ~ 3 ■ 10^® GeV. In minimal SUSY models, the values of 
all three gauge coupling constants gi coincide at this scale: gi{Mx) = gcvT 0- This 
relationship between the gauge coupling constants arises within Grand Unihed Theories 
0. In such theories, all observable gauge bosons and their superpartners (gaugino helds) 
belong to the same multiplet; therefore, the masses of all gauginos also coincide at the 
scale Mx- In the following, we everywhere set Mq,(Mx) = M 1 / 2 . The parameters M 1 / 2 , 
rri^, AafBj, and Bag dehned in this way at the Grand Unihcation scale should be treated 
as boundary conditions for the set of renormalization group equations that describes the 
evolution of the parameters of a soft breaking of supersymmetry down to the electroweak 
scale. 
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Within the minimal SUSY Standard Model (MSSM), an exact analytic solution to the 
renormalization group equations exists at tan ~ 1 [§] , in which case the Yukawa coupling 
constants hb and hr for the 6-quark and r-lepton, respectively, are negligibly small. 
This solution makes it possible to analyse the evolution of the f-quark Yukawa coupling 
constant ht{t) and of the parameters of a soft breaking of supersymmetry. For the 
coupling constant ht{t), this solution has the form 

m) 

where Yt{t) = h^/{A7rY and t = ln{Mx/Q'^). The explicit expressions for the functions 
Et{t) and Ft{t) are presented in the Appendix (see (|^ ). At low energies, the second term 
in the parentheses on the right-hand side of is much less than unity at sufficiently 
large values of hf{0); as a result, all solutions (H) to the renormalization group equations 
are focused in a narrow interval near the quasi-hxed point YQFp(f) = Et{t)/QFtit) 0. 
Formally, a solution of this type can be obtained by making Yf(0) tend to inhnity in (^. 

Along with the Yukawa coupling constant for the f-quark, solution to the renormalization 
group equations for the corresponding trilinear coupling constant At for scalar helds and 
the combination 971^ = mQ + m^ + m| of the scalar-particle masses approach the infrared 
quasi-hxed point with increasing h^(0). An analytic solution for these parameters can 
be represented as 



A (f\ — A I'n'i I \/T (tEt{t) Ft{t) 

A(t) -^ ^ 

<mut} = W(0) - .4?(0)) ll + (a,(0)|| - 

, X X V V 9 X N V O ^ V 




(5) 


X 




+ ^1/2 


± (.2E't{ty 

dt V Et{t) 


mit) 

Et{t) 


\ mJ 


where et{t) = Yt{t)/Yt{0). In the regime of strong Yukawa coupling, in which case 
ht{0) S> 5'gut( 0)) fhe natural small parameter et(t) arises in the theory being considered. 
In the infrared region, the dependence of the solutions in (^) and (|^) on the boundary 
conditions at the Grand Unihcation scale disappears almost completely. Near the infrared 
hxedpoint, we have Yt{to) ^ YQFp(fo), A{to) ~ ^qfp(^o), and^ 9JlQpp(fo), where 
AQFp(t) and 9?tQFp(^) are expressed in terms of only the gaugino mass at the scale Mx 
and to = 2 ln(Mx/Mf°^®), the pole t-quark mass being approximately equal to 

175 GeV. The deviations from YQFp(t), AQFp(t), and 911 qfp(^) are determined by the 
ratio €t{t)/Et{t) ^ 1/{6Ft{t)Yt{0)), which is of order l/(10/i^(0)) at the electroweak scale. 
The properties of solutions to the renormalization group equations within the MSSM and 
the spectrum of particles in the infrared hxed point regime at tan ~ 1 were investigated 

El- 


m 
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At large values of tan/3 (about 50 to 60), the Yukawa coupling constants hb and hr 
are of order ht] for this case, an exact analytic solution to the renormalization group 
equations has not yet been found. Nonetheless, the detailed investigations performed 
in 1,0 , ||14|| revealed that, in this region of the parameter space, solution to the set 
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of nonlinear differential equations being studied also approach the infrared quasi-hxed 
point, the basic properties of the solutions remaining unchanged. 


A reduction of the number of independent parameters in the vicinity of the infrared hxed 
point at tan/? ~ 1 make it possible to obtain, to a sufficiently high degree of precision, 
an upper limit on the mass of the lightest Higgs boson. In the case being considered, 
this upper limit does not exceed 94 ± 5 GeV 10.11 ,[|l5f|, the constraints on the mass of 


the lightest Higgs boson from LEP H [jT^ being such that a considerable part of solutions 
approaching the quasi-hxed point at tan /? ~ 1 have already been ruled out by the existing 
experimental data. This gives an additional incentive to study the Higgs sector and the 
renormalization group equations and their solutions within more involved SUSY models. 


Of extensions of the MSSM, the simplest one that makes it possible to preserve the 
unihcation of the gauge constants and which leads to a higher upper limit on the mass of 
the lightest Higgs boson is the nonminimal SUSY Standard Model (NMSSM) introduced 


m 


H.E 


The Higgs sector of this model contains an extra singlet superheld Y in 
addition to two doublets ffi and H 2 . The upper limit on the mass of the lightest Higgs 
boson within the NMSSM attains a maximum value in the regime of strong Yukawa 
coupling, in which case the Yukawa coupling constants U/O) are much greater than 
the gauge coupling constant dcuT = 9 g \ jt / parameter space region being 
considered, solutions to the renormalization group equations within the NMSSM are 
attracted to quasi-hxed lines or surfaces in the space of Yukawa coupling constants. In 
the limit Uj(0) —> 00 , all solutions to the set of differential equations in question are 
concentrated near quasi-hxed points |T^, which arise as the result of intersections of Hill 


lines or surfaces with the invariant line connecting the stable hxed point at Y, S> a* 
with the stable infrared hxed point within the NMSSM 
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For the parameters of a soft breaking of supersymmetry, the behaviour of solutions to 
the renormalization group equations within the NMSSM is studied here near infrared 
quasi-hxed points. Approximate solutions for the trilinear coupling constants Ai{t) and 
for the combinations 911^ (t) of the scalar particle masses are presented in the Appendix. 
In the regime of strong Yukawa coupling, where dcuT "C U/O) -C 1, Ai(t) and 9)T^(f) at 
the electroweak scale are concentrated near some straight lines or surfaces in the space 
of parameters being considered. With increasing Yi(0), the dependence of the trilinear 
coupling constants and of the combinations of the scalar particle masses on Aj(0) and 
becomes weaker; in the limit Yj(0) —>■ 00 , solutions to the renormalization group 
equations for the parameters of a soft breaking of supersymmetry approach a quasi-hxed 
point. 


2 Renormalization of the soft SUSY breaking terms 
for universal boundary conditions 

The superpotential in the NMSSM involves a large number of Yukawa coupling constants. 
At tan/? ~ 1, however, all of these are small, with the exception of the f-quark Yukawa 
coupling constant ht, the self-interaction constant x for the neutral scalar held Y, and 
the constant A characterising the interaction of the held Y with the doublets Hi and H 2 - 
In the regime of strong Yukawa coupling, the aforementioned constants can be chosen in 
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such a way that, at the scale Mx, the Yukawa coupling constant hb for the 6 -quark is 
equal to the Yukawa coupling constant for the r-lepton This relationship 

between hb and hr is usually realised in the minimal schemes of unification of gauge 
interactions 0 . 


Disregarding all Yukawa coupling constants, with the exception of h^, A, and x, we can 
represent the total superpotential within the NMSSM in the form 


X, 


W = + -Y^ + h,(H,Q)U‘^ 


( 6 ) 


By construction, the superpotential of the nonminimal SUSY model is invariant under 
the discrete transformations y'^ = of the Z 3 group IS. Second order terms in 


superhelds do not satisfy this condition; therefore, they have been eliminated from the 
superpotential ® . Upon a spontaneous breaking of gauge symmetry, the field Y develops 
a nonzero vacuum expectation value ((Y) = yl\/2). Thus, a mixing of the doublets Hi 
and H 2 , which is necessary for the emergence of the vacuum expectation value vi of the 
Hi doublet (without this vacuum expectation value, down quarks and charged leptons 
remain massless), is generated within the NMSSM. 

That the extra superfield Y, which is a singlet with respect to SU{2) ® U{1) gauge 
interactions, and the coupling constant A are introduced in the superpotential (||) leads 
to an increase in the upper limit on the mass of the lightest Higgs boson in relation to 
that within the MSSM: in the NMSSM, it is 135 GeV |^, which is 7 — 10 GeV greater 
than the corresponding value in the minimal SUSY model, the largest value of the mass 
of the lightest Higgs boson being attained in the regime of strong Yukawa coupling. 

Upon a soft breaking of supersymmetry due to SUGRA effects, scalar fields acquire masses 
m?; in addition, a trilinear constant Ai for the interaction of scalar fields is associated with 
each Yukawa coupling constant in the total Lagrangian of the theory. That the models 
being considered involve a large number of unknown parameters of a soft breaking of 
supersymmetry is one of the main flaws in such models. The hypothesis of universality of 
these constants at the scale Mx makes it possible to reduce their number in the NMSSM 
to three (the mass mo of all scalar particles, the trilinear coupling constant A for the 
interaction of scalar fields, and the gaugino mass M 1 / 2 ), whereby the analysis of the 
spectrum of SUSY particles is significantly simplified. Naturally, the universal parameters 
of supersymmetry breaking arise in the minimal SUGRA model (see 0,|^]) and in some 
string models @]. Even in the region of low energies, the hypothesis of universality 


of fundamental parameters makes it possible to avoid the emerging of flavour-changing 
neutral currents. Thus, the minimal set of the fundamental parameters in the NMSSM 
includes, in addition to the constants of the Standard Model, five unknown constants 
(A, X, A, mo, and Mi/ 2 )- Within the nonminimal SUSY model, the spectrum of the 
superpartners of observable particles and Higgs bosons for universal boundary conditions 
was studied in 


27 


The full set of renormalization group equations that describes the evolution of Yukawa 
coupling constants, the trilinear coupling constants Aj(t), and the scalar particle masses 
rnfit) within the nonminimal SUSY model from the Grand Unification scale down to the 
electroweak scale can be found in 10.11 


(see also Appendix). This set of equations is 


nonlinear even in the one-loop approximation; in view of this, its analytic solution has 
not yet been found. All equations that form the set being considered can be partitioned 
into two groups. The first includes equations that describe the evolution of gauge and 
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Yukawa coupling constants. If their evolution is known, the remaining equations from 
the set of renormalization group equations can be considered as a set of linear differential 
equations for the parameters of a soft breaking of supersymmetry. In solving this set of 
equations, it is necessary to integrate, hrst of all, the equations for the gaugino masses 
and for the trilinear coupling constants for the interactions of scalar particles. A general 
solution to the set of linear differential equations 

where the matrix Sij{t) and the column vector (nonhomogeneous term) Fi{t) are known, 
has the form 


t 

y,{t) = j (7) 

0 

where is a solution to the homogeneous equation d^ij(t)/dt = with 

the boundary conditions ‘hjj(O) = Sij. Since we have |/i(0) = (A, A, A) for Aj(f) if the 
fundamental parameters are chosen in a minimal way and since Fi{t) ~ M 1 / 2 , the trilinear 
coupling constants for the interactions of scalar helds are given by 


Aiit) = ei{t)A + fi{t)Mi/ 2 . (8) 

These solutions for Ai{t) must be substituted into the expressions on the right-hand 
sides of the renormalization group equations for the scalar particle masses, the functions 
Fi{t) involving terms proportional to A^, AM 1 / 2 , and Considering that, in the case 

being considered, ?/j(0) = (mg, tuq, mp, mg, rug), we can represent the required solution 
for mfit) as 


m-(f) = ai{t)ml + bi{t)Ml/^ + Ci{t)AMi /2 + di(t)A^. (9) 

The functions Ciit), fi(t), ai(t), bi(t), Ci(t), and di(t), which determine the evolution of 
Aj(f) and rn^it), remain unknown, since an analytic solution to the full set of renormaliza¬ 
tion group equations within the NMSSM is unavailable. These functions greatly depend 
on the choice of values for the Yukawa coupling constants at the Grand Unihcation scale 
Mx- At the electroweak scale t = to, relations (||) and (||) specify the parameters A^(fo) 
and m^(fo) of a soft breaking of supersymmetry as functions of their initial values at the 
Grand Unihcation scale. 

The results of our numerical analysis, which are presented in the Table, indicate that, 

with increasing Yi(0), where Yt{t) = Yx{t) = = 77^’ 

[Yk)^ (47r)"' (47r)"' 

functions ei(fo), q(^o), and diit ^) decrease and tend to zero in the limit Yi(0) — ^ cxd, 
relations d^) and (P) becoming much simpler in this limit. Instead of the squares of the 
scalar particle masses, it is convenient to consider their linear combinations 

^l{t) = ml{t)+ml{t)+ml{t), ( 10 ) 

ml{t) = 
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in analysing the set of renormalization group equations. In the case of universal boundary 
conditions, the solutions to the differential equations for DJtfit) can be represented in the 
same form as the solutions for rrifit) (see (D) ; that is 


+ Ci{t)AMi/2 + di{t)A^. 


( 11 ) 


Since the homogeneous equations for Ai{t) and (t) have the same form, the functions 
diit) and ej(f) coincide; in the limit of strong Yukawa coupling, the dependence 


disappears in the combinations (W) of the scalar particle masses as the solutions to the 
renormalization group equations for the Yukawa coupling constants approach quasi-fixed 
points. This behaviour of the solutions implies that Ai{t) and corresponding to 

hi(0) ^ dj(0) also approach quasi-fixed points. As was shown in |^, two quasi-fixed 
points of the renormalization group equations within the NMSSM are of greatest interest 
from the physical point of view. Of these, one corresponds to the boundary conditions 
1^(0) = Ta( 0) S> ttj(O) and Y^(0) = 0 for the Yukawa coupling constants. The hxed 
points calculated for the parameters of a soft breaking of supersymmetry by using these 
values of the Yukawa coupling constants are 


pX^ito) ^ 1.77, 
(^o) ^ -0.42, 


47 ( 

(to) ^ -2.^ 


p^/(fo) ~ 6.09, 

,QFP/ 


( 12 ) 


where pAiif) = Aj(t)/Mi /2 and Pm^{t) = Since the coupling constant x for 

the self-interaction of neutral scalar fields is small in the case being considered, A^{t) 
and do not approach the quasi-hxed point. Nonetheless, the spectrum of SUSY 

particles is virtually independent of the trilinear coupling constant A:^ since x —0. 


In just the same way, one can determine the position of the other quasi-fixed point for 
Aiit) and Tlfit), that which corresponds to Rxo = 3/4, = 3/8. The results are 


pX""(^o)^ 1.73, (to) ^-0.43, p:^^"(to) ~ 0.033, 
„QFP/'+w ^ no ^QFP,- 


^r( 

,QFPi 


(to) ~ 6.02, p^/ (to) ^-2.34, p^,^ (to) ^ 0.29, 


(13) 


where R\o = Ya(0)/17(0) and R^^ = Y>^(0)/Yt(0). It should be noted that, in the 
vicinities of quasi-hxed points, we have p^^(to) < 0. Negative values of 91t|(to) lead to 

a negative value of the parameter m^to) in the potential of interaction of Higgs helds 
(see Table). In other words, an elegant mechanism that is responsible for a radiative 
violation of SU{2) 0 U{1) symmetry and which does not require introducing tachyons 
in the spectrum of the theory from the outset survives in the regime of strong Yukawa 
coupling within the NMSSM. This mechanism of gauge symmetry breaking was hrst 


discussed in by considering the example of the minimal SUSY model. 


The evolution of the constants of a soft breaking of supersymmetry, Ai(t) and OJl^(t), is 
illustrated in Figs. 1-4 for Yi(0), Ya( 0) ^ d(0) and = 0, as well as for Y^(0) 3> d(0). 
Although the ratio A/M 1/2 is varied between —1 and 1, with the scalar particle mass mg 
lying in the range 0 < mg < solutions to the renormalization group equations are 

focused in a narrow interval at low energies. However, A^{t) and 91t^(t) are concentrated 
near zero, since the neutral field Y is not renormalized by gauge interactions; therefore, 
they remain dependent on the initial conditions (see also Table). The values of A*(to) and 
91Ij(fo) show the weakest dependence on A and mg, because the former are renormalized 
by strong interactions. 
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By using the fact that ^i(t) as determined for the case of universal boundary conditions 
is virtually independent of ml, we can predict ai{to) values near the quasi-hxed points. 
The results are 


1 ) RxO — 1 , R:M — 0 , 

1 4 5 

= O'ui.to) = ai{to) = ag{to) = 02 (^ 0 ) = ; 

2) R^o = 3/4, R^ = 3/8, 

^^(^0) 0 , ttl(to) ^2(^0) 2 ? ^q\to) g? ^uyto) — . 


(14) 


To do this, it was necessary to consider specihc combinations of the scalar particle masses, 
such as mlj—2mQ, mQ+m\j—ml+m\, and 2mf (at x = 0), that are not renormalized 
by Yukawa interactions. As a result, the dependence of the above combinations of the 
scalar particle masses on ml at the electroweak scale is identical to that at the Grand 
Unification scale. The data in the Table show that the predictions in (H) agree fairly 
well with the results of numerical calculations. 


3 Lines and planes in the space of parameters of the 
soft SUSY breaking 


By considering the example of an exact analytic solution Y;^ = 0, we can see how the A and 
ml dependence of the parameters of the parameters of a soft breaking of supersymmetry 
disappears. From this point, our analysis of solutions to renormalization group equations 
for Ai{t) and will not be restricted to the case of universal boundary conditions; 

that is, Aj(0) and fOf^(O) will be considered as independent boundary values. In the case 
being considered, the full set of renormalization group equations within the NMSSM can 
be partitioned into two subsets - of these, one coincides with the set of renormalization 
group equations within the MSSM, while the other describes the evolution of A^{t), 

and For Ya = 0, the evolution of At(t) and fOlfit) is determined by relations (H). 

For the other two parameters A,^{t) and of a soft breaking of supersymmetry, 

which correspond to the Yukawa coupling constant x, we obtain 

A^it) = A^{0)e^{t), 

mlit) = - Al{0)e^m - e.(t)), ^ ^ 


where e^(t) = Y^(f)/Y„,(0) = 1/(1 + 6Y^{0)t). With increasing Y><(0) and 4/(0), the 
values of the Yukawa coupling constants at the electroweak scale approach the quasi- 
fixed point, while solutions to the renormalization group equations for the parameters of 
a soft breaking of supersymmetry are attracted to the straight lines A^ = 1.67Mi/2 and 
A^ = 0 in the {At, A^) plane and to the straight lines = 5.49M/^2 ^-nd 911^ = 0 in the 
(91lj,9Jt^) plane. The dependence on the boundary values Aj(0) and 971^(0) disappears 
only in the limit Yj(0) —*• 00 , which corresponds to a quasi-hxed point in the {pt, p>t) 
plane. 


At nonzero values of Ya, it is possible to construct, on the basis of an approximate solution 
for the Yukawa coupling constants within the NMSSM |T^, approximate solutions to the 


renormalization group equations for the parameters of a soft breaking of supersymmetry. 






Near the quasi-fixed points given by (|^) and these solutions describe the evolution 
of At{t) and to a fairly high precision (about 1%). The relative deviations of the 

approximate solutions from numerical ones are significantly greater for A\{t) and 
At the electroweak scale, the relative error here is as large as 20 — 30% near quasi- 
fixed points. Finally, the approximate solutions for the parameters of a soft breaking 
of supersymmetry that correspond to the Yukawa coupling constant x are characterised 
by the poorest accuracy, providing only a qualitative description of the behaviour of the 
numerical solutions. Our results lead to the conclusion that the relative deviation of the 
approximate solution being studied from the precise numerical solution decreases with 
increasing contribution of gauge interactions to the renormalization of the corresponding 
parameter of a soft breaking of supersymmetry. 


Approximate solutions to the renormalization group equations within the NMSSM for 
the trilinear coupling constants and the combination ([l^ of the scalar particle masses 
are presented in the Appendix. In the regime of strong Yukawa coupling at Aj(0) = A 
and nJt^(O) = 3mg, the dependence of these solutions on the initial conditions vanishes in 
proportion to 1/Yj(0). Disregarding (9(1/1^(0)) terms, we find, for nonuniversal boundary 
conditions, that 


/ At{t) \ RxoFxjt) f AtiO) - Aa(0) \ 

+ 2RxoFx{t) VnjiKo) - +''' ’ 

3Ft{t) / Aa(0)- A(0) \ 

^6R{t) + 2RxoFx{t) Vmt2(o) - mm) 

2R^t / Aa(0 )-A^(0 ) \ ^ ’ 

6RxoFx{t) + 6R^t Wi(0) - mm) ’ 
fAm \ QRxoFxit) f Am - ^a(O) \ 

[mat]) ^QR^xit) + 6R^t Vafi?,( 0 ) - mm) 

In (|l6l), we discarded terms proportional to M 1 / 2 , Ai{0)Mi/2, and Aj(0)Aj(0). The 

explicit expressions for the functions Fi{t) are presented in the Appendix (see (^)). 
From the results that we obtained, it follows that, in the case of nonuniversal boundary 
conditions, the approximate solutions to the renormalization group equations for Ai{t) 
and m‘f{t) remain dependent on the difference of the values of these constants at the 
Grand Unification scale. At the same time, the linear combinations 


9Ft{t)At{t) + 3RxoFx(t)Ax(t) + R^tA^ = const, 
9Ft{t)ml{t) + 3RxQFx{t)m\{t) + R^tmm) = const', 


of the parameters undergo virtually no changes in response to variations of Aj(0) in the 
first combination and of m‘f{0) in the second combination. At the electroweak scale 
(t = to), relations (p!?]) specify straight lines (at = 0) and planes (at h%(0) S> d(0)) in 
the space of the parameters of a soft breaking of supersymmetry. In the regime of strong 
Yukawa coupling, the approximate solutions to the renormalization group equations are 
attracted to these straight lines and planes. 

The results of our numerical analysis, which are illustrated in Figs. 5 and 6, indicate 
that, in the vicinity of the infrared fixed point at = 0, solutions to the renormalization 
group equations at the electroweak scale are indeed concentrated near some straight lines 
for the case where the simulation was performed by using boundary conditions uniformly 
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distributed in the {At, A\) and the OJl^) plane. The strength with which these 
solutions are attracted to them grows with increasing Ti(0). The equations for the lines 
being considered can be obtained by htting the numerical results displayed in Figs. 5 
and 6. This yields 

+ 0.147(0.121)Aa = I. 7 OM 1 / 2 , g 

+ 0.147(0.121)971^ = ^ ^ 

For Y^{0) S> di(0) solutions to the renormalization group equations are grouped near 
planes in the space of the parameters of a soft breaking of supersymmetry {At, Ax, A,^) 
and {TtlTtlml) (see Figs. 7-9): 


At + 0.128(0.091)71a + 0.022(0.0105)71^ = I. 68 M 1 / 2 , 
+ 0.128(0.091)971^ + 0.022(0.0105)971^ = 5.77Ml/^. 


In (1^) and (pT[) , the predictions for the coefficients of Ax, A^, 97l|, and 971^ according 
to the calculations at Rxo = 1, R>to = 0 and at Rxo = 3/4, = 3/8 on the basis of 

the approximate solutions (^) and (p!7D to the renormalization group equations within 
the NMSSM in the regime of strong Yukawa coupling are given in parentheses. It can 
be seen from Figs. 7 and 8 that, as the values of the Yukawa coupling constants at the 
Grand Unification scale are increased, the areas of the surfaces near which the solutions 
Ai{t) and 97t^(t) are concentrated shrink in one of the directions, with the result that, at 
Yi(0) ~ 1, the solutions to the renormalization group equations are attracted to one of 
the straight lines belonging to these surfaces. 


Thus, the approximate solutions presented in this study lead to qualitatively correct re¬ 
sults. However, an analysis of numerical solutions to the renormalization group equations 
revealed that, with increasing 1/(0), only in the regime of infrared quasi-fixed points (that 
is, at Rxo = 1, R^ = 0 or at Rxo = 3/4, R,^ = 3/8) do ei{to) and di{to) decrease in 
proportion to 1/Yj(0). Otherwise, the dependence on A and tUq disappears much more 
slowly with increasing values of the Yukawa coupling constants at the Grand Unification 
scale - specifically, in proportion to (^(0))“*^, where 5 < 1 (for example, 6 = 0.35 — 0.40 
at X = 0). In the case of nonuniversal boundary conditions, only when solutions to 
the renormalization group equations approach quasi-fixed points are these solutions at¬ 
tracted to the fixed lines and surfaces in the space of the parameters of a soft breaking 
of supersymmetry; in the limit 1/(0) —>• 00 , the parameters Ai{t) and 97t^(t) cease to be 
dependent on the boundary conditions, in contradiction with the results presented in 16 
(|l^). Finally, we note that, by using approximate solutions to the renormalization group 
equations within the NMSSM, it is impossible to explain the emergence of those lines in 
the {At, Ax, A,^) and (971^, 971^, 971^) spaces to which solutions of the equations in question 
are attracted when all three Yukawa coupling constants are 1/(0) ~ 1. 


4 Analysis of solutions of the renormalization group 
equations near quasi—fixed points 

That, in the regime of strong Yukawa coupling, quasi-fixed points, lines, and surfaces 
appear in the space spanned by the parameters of a soft a breaking of supersymmetry is 
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explained by fast changes in the Ynkawa conpling constants at the initial stage. In view 
of this, the entire scale from Mx to can be partitioned into two nneqnal intervals. 

Within the initial stage corresponding to Mx > O' > O'! ~ 10^^ GeV, the Ynkawa conpling 
constants are signihcantly greater than the gange conpling constants and decrease fast 
with increasing t. In this region, the evolntion of all fnndamental parameters of a soft 
breaking of snpersymmetry is determined by Yi(t), so that all gange conpling constants 
ai{t) can be disregarded for a first approximation. At the initial stage of evolntion, the 
one-loop renormalization gronp eqnations for Yi{t) can then be represented in the form 

^c„Yi(t)y ( 20 ) 


where Cij > 0. If the set of (|^) has a fixed point, all the Ynkawa conpling constants are 


proportional to each other in the vicinity of this point. In analysing o. it is therefore 
convenient to introdnce the additional constant Yq satisfying the eqnation 


dYojt) 

dt 




and to investigate, instead of Yi(t), the ratio rj(t) = Yi{t)/Y^it). The fixed point of 


the renormalization gronp eqnations (1^) is then determined by solving the set of linear 
algebraic eqnations 


dij) — 0 . ( 21 ) 


Linearising the set of renormalization group (|^) near the fixed point ri{t) = + 

we obtain 


dOi 

dz 


dAi 

dz 

dmj 

dz 


j 

Y CijVjOjAj'j , 
j 

^ (cpT° [m] + Afj + Cijrpj [m] + Aj)) , 


( 22 ) 


where 2 ; = ln(Yo(^)/ho(0)) = — ln(l + Y^t). A general solution to this set of linear 
differential equations at = 0 has the form 


M^) = '^akUikcY^, 

} (23) 

= '^PkUikeY^ + '^XkeYMikU^^ j Al{y)e 

k k,n Q 

where Xk are the eigenvalues of the matrix CijVj, Uik is the eigenvector corresponding to 
the eigenvalue A^, and ak and (3k are constants of integration that are expressed in terms 
of Aj(0) and 071^(0). It should be noted that all eigenvalues \k are positive. They are 
determined by solving the characteristic equation det(cjjr° — X'Sij) = 0. 
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From (1^) , it follows that, at A' = 1, one of the columns of the matrix {cijr^ — X'Sij) 


is a linear combination of all other columns; hence, its determinant is identically equal 
to zero. This means, that the value of Aq = 1 is a root of the characteristic equation. 
The eigenvector Uio = (1,..., 1) corresponds to this eigenvalue. As a rule, all other roots 
of the characteristic equation are less than unity. With increasing t, the dependence of 


solutions (1^) on the boundary conditions at the Grand Unihcation scale becomes weaker. 
At g ~ 10^^ — 10^^ GeV, in which case the Yukawa coupling constants are on the same 
order of magnitude as the gauge coupling constants, the values of Ai{z) and at 


this scale that are calculated by formulas (P3| ) should be treated as new boundary values 
for the full set of renormalization group equations within the NMSSM, which is presented 
in the Appendix. 

In the case where only the f-quark Yukawa coupling constant is nonzero, the set of 
differential (1^) for the Yukawa coupling constants reduces to one equation; as a result, 
the full nonlinearised set of renormalization group equations can be solved exactly in the 
limit dj —> 0. The results are given by 


y.it) = 


y.(o) 


Atit) — At{0)et{t), 


l + 6Yt(0)G 

- A?(0)6,(f)(i - 


(24) 


where et{t) = Yt(f)/Yt(0). The solution in (^41) coincides in form with solution (p3|) to 
the above linearised set of differential equations. Formally, it can be obtained by setting 
Afi /2 = 0 and E{t) = 1 in the exact solution to the NMSSM renormalization group 
equations that is given by (§) and (§). A comparison of the analytic solution (^41) in the 
regime of strong Yukawa coupling with the exact solution to the renormalization group 
equations within the minimal SUSY model reveals that, as the result of the evolution of 
At{t) and from the scale qi ~ 10^^ GeV to the electroweak scale, there arises, in 

addition to et(t), an extra factor 1/E(t) in front of At(0) and 9It4(0) in (^. This new 
factor is independent of the initial conditions for the f-quark Yukawa coupling constant 
at the Grand Unihcation scale. 


From relations (p3D, it follows that an expansion of solutions to the renormalization group 
equations in powers of the small parameters et{t) - this expansion emerges in the regime 
of strong Yukawa coupling within the MSSM - is also valid in more involved models. 
At the electroweak scale, solution (^) to the renormalization group equations can be 
represented as 


Mt) \ 

mm) 




(25) 


where the quantities Vikit) are weakly dependent on the Yukawa coupling constants at 
the Grand Unihcation scale and et{t) = Yt(f)/V(0). In expressions we have dis¬ 
carded terms proportional to M 1 / 2 , and A^. The extra factors Vikit) 

appear upon renormalizing the parameters of a soft breaking of supersymmetry from 
q ~ 10^^ — 10^^ GeV to g ~ m^. At x = 0, there exists, in addition to Aq = 1 
and Uio = (1,1), yet another eigenvalue Ai = 3/7 and the corresponding eigenvec¬ 
tor Uii = (1,-3), whose components specify {At, Ax) and {DJlt,VJl\). With increasing 
17 ( 0 ) = Y\(0), the dependence on cxq and Pq becomes weaker and the solutions at f = to 
are concentrated near the straight lines {At{ai), Aa(q;i)) and (0Jl^(/3i), DJl\{Pi)). In order 
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to obtain the equations for these straight lines, it is necessary to set Aa( 0) = —3y4t(0) and 
= —3971^(0) at the Grand Unihcation scale. At the electroweak scale, there then 
arise a relation between At{to) and A\{tQ) and a relation between Tlfito) and 


At + 0.137Aa = I. 7 OM 1 / 2 , 
+ 0.137971^ = 


These relations agree well with the results obtained above and presented in (^). More¬ 
over, the equations deduced for the straight lines at 17 ( 0 ) ~ 1 by htting the results of the 
numerical of the numerical calculations nearly coincide with those in (l2^). 


When the Yukawa coupling constant x is nonzero, the matrix has three eigen¬ 


values. One of these is Aq = 1, while the other two are Ai = 


3 + 75 
9 


0.58 and 


3 — 75 

A 2 = —^~ 0.085. The corresponding eigenvectors Un and Ui 2 , which specify 
{At,Ax,A^) and (971?, 97l|, 97t^) are given by 




a(l + V5)' 

+5 


Ui2 = 



(27) 


An increase in Ta(0) = 2Y^(0) = ^^(0) leads to the following: first, the dependence of 

Ai(t) and 97t?(f) on oq and Pq disappears, which leads to the emergence of planes in the 
space spanned by the parameters of a soft breaking of supersymmetry. 


At + 0.103Aa + 0.0124A^ = 1.69Mi/2, 
971? + 0.103971? + 0.0124971?, = 5.78M?/2; 


(28) 


after that, the dependence on ai and Pi becomes weaker at 17(0) ~ 1. This means 
that, with increasing initial values of the Yukawa coupling constants, solutions to the 
renormalization group equations are grouped near some straight lines and we can indeed 
see precisely this pattern in Figs. 7-9. All equations presented here for the straight lines 
and planes in the 971? space were obtained at Aj(0) = 0. 


From relations (^6|) and (p8|) , it follows that At(to) and 97l?(fo) are virtually independent 
of the initial conditions; that is, the straight lines and planes and orthogonal to the At 
and 971? axes. On the other hand, the A^(to) and 97t?,(to) values that correspond to 
the Yukawa self-interaction constant Y for the neutral fields are fully determined by the 
boundary conditions for the parameters of a soft breaking of supersymmetry. For this 
reason, the planes in the {At, Ax, A,P) and (971?, 971?, 971?,) spaces are virtually parallel to 
the Aye and 971?, axes. 


In the case of nonuniversal boundary conditions, the parameters of a soft breaking of 
supersymmetry as functions of Aj(0) and 97t?(0) in he regime of strong Yukawa coupling 
are determined by (ei(f))^™ (see (^)), where Amin is the smallest eigenvalue of the matrix 
CijV^ and epto) ~ l/h?(0). If the values of all Aj(0) and all 97t?(0) coincide, a solution to 
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be linearised set of differential equations at 7 ^ 0 can be represented in the form 




A,{z) =A(0) + (e^ - 1) 


dOiiz] 


ml{z) =im 2 ( 0 ) e" + (e" - 1 ) 


dz 
dOAz) 


dz 


(29) 


+ 2170) (e"-l) 


x(e' + l(e' + l)^5^ + (e'-l)^ 

\ 2 dz dz^ 

where 7 ^ are constants of integration. In the stable hxed point regime, where all 9i are 
equal to zero ,the evolution of Ai{t) and is identical to that within the MSSM 

(see (1^). In the limit being considered, the dependence on ^4(0) and II)T^(O) dies out 
fast, in proportion to Ytif )and the parameters Ai{t) and Tlj{t) approach the 
hxed point. Since 6i ~ > — cx), the dependence of the parameters of a 

soft breaking of supersymmetry on the initial conditions is determined by in 

all other cases - that is, by the strength with which the quasi-hxed point attracts the 
solutions to the renormalization group equations for the Yukawa coupling constants in 
the limiting case of I^j(O) di(0). As a rule, Amin is signihcantly less than unity; for 
example, Amin = 3/7 ~ 0.43 at x = 0 and Amin ~ 0.085 at nonzero x. This is the reason 
why, in the regime of strong Yukawa coupling, the dependence on A and mg dies out at 
a rate much lower than that which is predicted by approximate solutions. Nonetheless, 
expressions (p^)-(|25|) and ( p9[) for Ai{t) and tend to zero in the limit 4/(0) —>• 00 ; 

for this reason, the Aj( 0 ) and 971^(0) dependence of the solutions being studied disappears 
completely. 


5 Conclusion 

An increase in the Yukawa coupling constants at the Grand Unihcation scale results in 
that, for Yt{t), Yx{t), and Y^{t), the solutions to the renormalization group equations 
within the NMSSM are attracted to infrared quasi-hxed points. We have shown that, 
as Yi{t) approach these points, the corresponding solutions for the trilinear coupling 
constants Ai{t) characterising scalar helds and for the combinations 911^ (t) of the scalar 
particle masses (see O) cease to be dependent on their initial values at the scale Mx 
and, in the limit 1 /( 0 ) —> 00 , also approach the hxed points in the space spanned by the 
parameters of a soft breaking of supersymmetry. Since the set of diherential equations 
for Ai{t) and rnf{t) is linear, the A, M 1 / 2 , and mg dependence of the parameters of a 
soft breaking of supersymmetry at the electroweak scale can be explicitly obtained for 
universal boundary conditions. It turns out that, near the quasi-hxed points, all Ai{t) 
and all 971^(t) are proportional to M 1/2 and respectively. Thus, we have shown 

that, in the parameter space region considered here, the solutions to the renormalization 
group equations for the trilinear coupling constants and for some combinations of the 
scalar particle masses are focused in a narrow interval within the infrared region. Since 
the neutral scalar held Y is not renormalized by gauge interactions, Aj^{t) and 9It^(t) are 
concentrated near zero; therefore they are still dependent on the initial conditions. The 
parameters At(tg) and DJl^ito) show the weakest dependence on A and mg because these 
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parameters are renormalized by strong interactions. By considering that the quantities 
dJUfito) are virtually independent of the boundary conditions, we have calculated, near 
the quasi-hxed points, the values of the scalar particle masses at the electroweak scale. 

In the general case of nonuniversal boundary conditions, the solutions to the renormal¬ 
ization group equations within the NMSSM for Ai(t) and are grouped near some 

straight lines and planes in the space spanned by the parameters of a soft breaking of 
supersymmetry. Moving along these lines and surfaces as Bi(0) is increased, the trilinear 
coupling constants and the above combinations of he scalar particle masses approach 
quasi-hxed points. The formation of the straight lines and planes can be traced by con¬ 
sidering the examples of an exact (at Yx = 0) and an approximate (for Yx ^ 0) solution. 
Although the approximate solution for the parameters of a soft breaking of supersymme¬ 
try in a nonminimal SUSY model (this solution has been obtained here for the hrst time) 
describes the evolution of A^{t) and very poorly, it leads to qualitatively correct 

results. 

However, our analysis of the behaviour of solutions to the renormalization group equations 
within the NMSSM in the regime of strong Yukawa coupling has revealed that the Aj(0) 
and (0) dependence of the trilinear coupling constants and the above combinations 
of the scalar particle masses dies out quite slowly, in proportion to (ei(t))^“‘“, where 
etit) = y)(t)/Yi(0) and Amin is the minimal eigenvalue of the matrix Cj^r^. As a rule, 
Amin is positive and much less than unity. For example, Amin = 3/7 at = 0 and 
Amin ~ 0.0085 at ^ 0. The above is invalid only for the solutions Ai{t) and ^f(t) 
that correspond to universal boundary conditions for the parameters of a soft breaking of 
supersymmetry and to the initial values of Rxo = 1, R>co = 0 and Rxo = 3/4, R^ = 3/8 
for the Yukawa coupling constants at the Grand Unihcation scale. They correspond 
to quasi-hxed points of the renormalization group equations for Yj(t). As the Yukawa 
coupling constants are increased, such solutions are attracted to infrared quasi-hxed 
points in proportion to et(t). 

It is precisely because of the inequality Amin 1 that straight lines in the (A^, Ax, A^ 
and (971^, !DT|, spaces, rather than infrared quasi-hxed points, play a key role 
in the analysis of the behaviour of solutions for Aj(t) and ^f{t) in the case where 
Yt(0),YA(0),Y.(0) ~ 1. In the space spanned by the parameters of a soft breaking of 
supersymmetry, these straight lines he in the planes near which Ai{t) and 971^ (t) are 
grouped in the regime of strong Yukawa coupling at the electroweak scale. A method has 
been proposed in the present study for deriving equations that describe the straight lines 
and planes being discussed. The straight lines and planes that were obtained by using 
this method or by htting the results of our numerical calculations are nearly orthogonal 
to the At and 971^ axes. This is because the constants At(to) and 971^(to) are virtually 
independent of the initial conditions at the scale Mx- On the other hand, the parameters 
A,^(to) and 97l^(to) are determined, to a considerable extent, by the boundary conditions 
at the scale Mx- At Rxo = 3/4 and R,^ = 3/8, the planes in the (At, Ax, A,^) and 
(ml ml, ml) spaces are therefore parallel to the A^ and 971^ axes. 

Within the NMSSM involving a minimal set of fundamental parameters, we have been 
unfortunately unable to construct, at U(0) 3> Q;j(0), a self-consistent solution that, on 
one hand, would lead to a heavy spectrum of SUSY particles and which, on the other 
hand, would make it possible to obtain, for the lightest Higgs boson, a mass value greater 
than that in the MSSM. In order to obtain such a solution, it is necessary to modify the 
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nonminimal SUSY model. In the near future, we plan to investigate the spectrum of the 
superpartners of observable particles and Higgs bosons within the simplest modihcation 
of the NMSSM such that it leads to a self-consistent solution in the regime of strong 
Yukawa coupling. 
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Appendix 


Set of the renormalization group equations of the NMSSM and an approxi¬ 
mate solution of it. 


The full set of one-loop renormalization group equations of the NMSSM can be repre¬ 
sented in the form 


dai 

dt 

dt 

dt 

d^ 

dt 

dMi 

dt 

dAt 

dt 

dt 

dA^ 

dt 

dm\ 

dt 

dm^ 

dt 


drriy 

dt 


16 13 

~ Yti^x + GYt —— 3(52 ~ ~(5i), 

3 15 

— Yx(4Yx + 2Y^ + 3Yt — 3(52 ~ 7 ^ 1 ), 

5 

-6Y^{Yx + Y^), 

— bidiiMi, 

16 13 

~ Y\Ax — QYtAt + —dsM^ + 3(52^2 -|- —diMi, 

3 15 

— AYxAx — 2Yj^A:^ — 3YtAt 302^^2 + ■rdiiffi, 

5 

— 61a-4a - QY^A^, 

3 

— Yxirn^ -|- 7T7.2 -f rn.^ -f Af^ -|- 35t2M2 + zrdiAf^, 

^ 5 

— Yx{m\ + ml + ml + Al) - 3Yt{ml + mfj + tuq + A'^) -f 

3 

-|- 3d'2M2 —diMi, 

5 

— 2Yx{ml + ml + ml + Al) - 2F^(3mJ -f A^), 


(30) 


dmfj 

dt 

dml 

dt 


1 fi 1 fi 

= - + mlj + ml + A^) + -d,Ml + -d,Ml 

16 1 

= — Yfi^ml mfj 7T7.Q -|- AI) H— —d^M^ Y 302-^2 Y -—diMl, 

O -LO 


where the index i runs through the values from 1 to 3, bi = 33/5, 62 = 1, ^3 = ~3, 
di{t) = ai{t)/{47r) = gf(t)/(Air)^, Yt{t) = hl{t)/{47r)‘^, Yx{t) = X‘^{t)/{47r)^, and 
= x^(f)/( 47 r)^. The variable t is defined in a conventional way: t = ln(M|'/Q^), 
where Mx = 3 ■ 10^® GeV. On the right-hand sides of (^), we have discarded terms 


proportional to the Yukawa coupling constants Yf, and Y^-, because their contribution is 
negligibly small at tan/? ~ 1. For the Yukawa and gauge coupling constants, the set of 
two-loop renormalization group equations within the NMSSM is presented in [^. Using 
the results obtained from an analysis of renormalization of various types of constants 
in SUSY theories and also the two-loop equations for the Yukawa and the gauge 


coupling constants, one can easily obtain the relevant renormalization group equations 
describing the evolution of parameters of a soft breaking of supersymmetry. Moreover, 
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the full set of two-loop renormalization group equations of the NMSSM can be obtained 
directly from the results presented in [^|. 


Although an exact analytic solution to the set of one-loop renormalization group equa¬ 
tions (pOD of the NMSSM is not known at present, an approximate solution for the 
parameters of a soft breaking of supersymmetry can be constructed on the basis of the 


approximate solution for the Yukawa coupling constants jlll : 


Ytit) 

Yxit) 

Y^{t) 

Etit) 

Exit) 

aiit) 


YtiO)E,it) 

(1 + 6Yi(0)F,(f))V2(i + QYtiO)Ftit) + 2Y,iO)F^it))Y^ ’ 

YxiO)Fxit) 

(1 + 6YtiO)Ftit) + 2YxiO)Fxit))Y^il + 6YxiO)Fxit))y^ 

1 

"" (l + 6YA(0)FA(f)+6Y^(0)t)V3’ 

__ 

(l + 6YA(0)FA(f) + 6n.(0)f)’ 


asity 

16/9 

'a2ity 

-3 

'aiity 

d(0) _ 


_d(0)_ 


Q!(0) 


-13/99 


, F^it) = / Etit’)dt', 


02 (t) 


-3 


aiit) 


a(0) a(0) 

d(0) 


-1/11 


, Exit) = / Exit’)dt’, 


1 -|- 6jd(0)t 


(31) 


To do this, it is necessary to replace the quantities cij and Yj in the approximate solution 
(0) by d' and Y/: 

a[ = ai{l + MiTj + Mif] 2MiMirjff) , 

Y/ = Yi{l + AiT] + Aifj + 2AiAir]fj + , (32) 

^ V = 


where 6 are Grassmann variables, and to construct an expansion in rj. This method for 
deducing solutions made it possible to simplify considerably the expressions for Atit) 
and 2iytfit) within the MSSM at tan/? r\j 1 and to obtain and analyse approximate 
solutions for the parameters of a soft SUSY breaking within the MSSM for tan/? S> 1 
33. In the most general case of nonuniversal boundary conditions, we have 


^iit) — Cjj(t)Aj(0) -|- fiit)Mi/2, 
j 

(^) ~ (®q(^)^/(0) + Cjj(f)Aj(0)Mi/2) + dijkit)AjiO)AkiO) + biit)My2y 

j jk 


where Cijit) = aijit), since the homogeneous equations for Aj(f) and OJl^(f) coincide, and 
where the functions e^/f), dijit), /j(f), hit), Cijit), and dykit) are expressed in terms of 
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the partial derivatives of the Yukawa coupling constants as 


Cijit) = Yj(0) 


d\nYi{t) 
dY,{0) ’ 


fi{t) = q;(0) 


bi{t) = 


d 


d\nYi{t) 

9d(0) 


9d(0) 

Cij{t) = 2d(0) 


d'(0) 
d 


da{0) 


d\nYi{t) 
da{0) 

dlnYiit) 


Y,(0) 


dYM 


= fi{t) + a(0) 
= 2d(0) 


dijk(t) ^^(0) 


d 


dYkiO) 


Y,(0) 


d\nYi{t) 

dY,{0) 


dfijt) 

5(5 (0) ’ 

deij{t) 

da{0) ’ 

= n(o) 


deij{t) 

dY,{oy 


( 34 ) 


Substituting the above approximate solutions for the Yukawa coupling constants within 
the NMSSM into (0), we obtain: 


e\\(t) 

(^) 

ext{t) 

e^xit) 

Cprf(t) 

etx{t) 

euit) 

hit) 

fht) 

Mt) 


^ 2YxiO)Fxit) _ YxiO)Fxit) 

1 + 6Ya(0)-Fa(^) + 6Y^(0)t 1 + 6 Ya(0)-Fa(^) 

_ YxiO)Fxit) 

l + 6Yi(0)F,(f) + 2YA(0)FA(t)’ 

_ 2Y^{0)t 

l + 6YA(0)FA(t) + 6n.(0)t’ 

_ miO)Ftit) 

l + QYt{0)Ft{t) + 2Yxi0)Fxit)^ 

6Yx{0)Fx{t) 

l + 6Yxi0)Fxit) + 6Y^{0)F 
6Y^{0)t 

l + 6Y,{0)F,{t) + QY^{0)F 

eUt) = 0, (35) 

_bA(0)FA(t)_ 

l + 6Yt(0)Fi(t) + 2YA(0)FA(t)’ 

3Yi(0)Fi(f) 3Yt(0)Fi(t) 

1 + 6Yi(0)Fi(f) 1 + 6Yi(0)Fi(t) + 2YA(0)FA(t) ’ 

tEht) YximExit) - Exit)) 2YxmtExit) - Exit)) 

E,{t) 1 + 6YA(0)FA(f) 1 + 6YA(0)FA(t) + 6Y^(0)f 

3Yt{0){tEt{t) - Ftjt)) + YA(0)(t^A(t) - Exit)) 
l + 6Yt{0)Ft{t)+2Yxi0)Fxit) 

6Yx{0)itE^it) - F,{t)) 
l + QYxiO)Fxit) + 6Y^iO)E 
tE'{t) 3YyO){tEtit) - Fyt)) 

Eyt) l + 6Yi(0)Ft(f) 

3YtmtEtit) - Ft{t)) + YA(0)(t^A(t) - FA(t)) _ 
l + 6Yi(0)Ft(f)+2YA(0)FA(f) 
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bx{t) 


b^{t) 

hit) 


cxxit) 


cx>cit) 

cxtit) 


c^xit) 

^3-C3<it) 

Cjtt it) 

hxit) 


cuit) 


± , GY^itE^jt) - F,it)f 

dt\ Exit)) l + QYxiO)Exit) il + 6YxiO)Exit))^ 

3YtiO)t^E',{t) + 2Y^iO)t^E'^it) 

1 + 6YtiO)Etit) + 2n(0)FA(t) 1 + 6n(0)FA(t) + 6Y^iO)t 

2i3YtiO)itEtit) - Etit)) + YxiO)itExit) - Exit)))^ 

(l + 6Ft(0)F,(t) + 2n(0)F;,(t))2 
UY^itExit) - Exjt))^ 
il + 6YxiO)Exit) + 6Y^iO)t)^' 

QY,iO)t^E',it) 36Y^iO)itE,it) - E^jt))^ 

1 + 6n(0)F,(t) + 6Y^iO)t (1 + 6YxiO)Exit) + QY^iO)t)^ ’ 

d_ (2Em\ _ mo)t^m) , isYmmit) -E,it))^ 

dt\ Etit)) l + QYtiO)Etit) il + 6YtiO)Etit))^ 

2i3YtiO)itEtit) - Etjt)) + YxmtExjt) - Exjt)))'^ 
il + 6YtiO)Etit) + 2YxiO)Exit))^ 

3YtiO)t^Elit) + YxiO)t^E'^it) 
l + 6YtiO)Etit) + 2YxiO)E,ity 

AYxiO)itExit) - Exit))il + 6Y^iO)t) 2YxiO)itExit) - Exit)) 
(1 + QYyO)Eyt) + 6F^(0)t)2 (1 + 6FA(0)F,(t))2 

YtiO)YxiO)Exit)itEtit) - Etjt)) 
il + QYtiO)Etit) + 2YxiO)Exit))^ 

Y,iO)itE,it) - E,it))il + 6YtiO)Etit)) 
il + QYtiO)Etit) + 2YxiO)Exit))^ ’ 

2AY^iO)YxiO)itEyt) - Exit))t 
(l + 6FA(0)F,(t) + 6F^(0)t)2 ’ 
^^ YtiO)YxiO)Etit)itExit)-Exit)) 

il + QYtiO)Etit) + 2YxiO)Exit))^ 

YtiO)itEtit) - Etit))il + 2n(0)FA(t)) 
il + QYtiO)Etit) + 2YxiO)Exit))^ ’ 

Y,iO)itE,it) - E,it))il + 6Y^m 
il + QYyO)Eyt) + QY^iO)t)^ ’ 

72Y^iO)YxiO)itExit) - Exit))t 

(l + 6FA(0)F,(t) + 6F^(0)t)2 ’ 

Ctj^it) 0, 

YtmxmxjtEtjt) - Ftjt)) 

il + QYtiO)Etit) + 2YyO)Eyt))^ 

YxmtExjt) - Eximi + 6YtiO)Etit)) 

il + 6YtiO)Etit) + 2YmFxim ’ 

6YtiO)itEtit) - Etjt)) ^ mtiO)YxiO)Etit)itExit) - Exit)) 

(1 + 6YtiO)Etim (1 + QYtimtit) + 2YxiO)Exim 
YmtEtjt) - Etit))il + 2YmFxit)) 

il + QYtiO)Etit) + 2YxiO)Exim ’ 
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dxxt{t) 

dx>c}i(i^ 

dxtt{t) 

d^xx{t) 

d}c}i}c(t') 

dtxx{t) 

dtxt{t) 

dttt (^) 
dx>d{t) 

djitt (^) 


2Yx{0)Fx{t){l + 6Y^{0)t) Yx{0)Fx{t) 

(1 + 6yA(0)FA(t) + QY^{0)ty (1 + 6Yx{0)Fx{t)y 
Yxi0)Fxmi + 6Yt{0)Ftit)) 
il + QYtiO)Ft{t) + 2YxiO)Fxit)y^ 

T ... ^ UY^mmtFxjt) 

A-Al J {i + QYx{0)Fx{t) + QY^{0)ty^ 

7 6Y,{0)Yx{0)F,it)Fx{t) 

AtAl J {l + QYt{0)Ft{t) + 2Yx{0)Fx{t)y' 

2Y^{0)t{l + 6Yx{0)Fx{t)) 
il + QYxi0)Fxit) + 6Y^{0)ty^ 

3Yt{0)Ftit){l + 2Yx{0)Fx{t)) 
il + QYtiO)Ftit) + 2YxiO)Fxit))^^ 

6Yx{0)Fx{t)il + 6Y^{0)t) 
il + QYxiO)Fxit) + QY^{0)ty' 

~ _ 36lA(0)y>r(0)-FA(t)t 

) - (l + 6FA(0)FA(t) + 6n.(0)t)2’ 

6n.(0)t(l + 6yA(0)i"A(t)) 

{l + QYx{Q)Fx{t) + YMt?' 

Yx{Q)Fx{t){l + QYmFt{t)) 

{l + QYmFt{t) + 2Yx{Q)Fx{t)r 

7 6y,(0)yA(0)F^(t)FA(t) 

^ (l + 6Ft(0)Fi(t)+2FA(0)FA(t))2’ 

3y,(0)F,(t) 3Y,{Q)Ft{t){l + 2yA(0)i"A(t)) 

(1 + QYmFt{t)Y (1 + QY^Ftit) + 2FA(0)FA(t))2 ’ 

d^\{t^ d^y^{t^ 0, 

dtx^tit^ dt:K\{fy dij^;^{t^ dt^{t) dtt^{t^ 0* 


(H) 


At Va( 0) = 0, the solution given by (|^) and ([^ coincides with the exact analytic solution 
of the set of renormalization group equations (^). Summing the above expressions for 


and 


ei{t) = ai{t) = Ci{t) = 


di{t^ ^ djjk (^) ? 

jk 


we find the functions ei{t), Ci{t), and di{t), which determine the dependence of the param¬ 
eters of a soft SUSY breaking on A, M 1 / 2 , and tuq for the universal boundary conditions. 
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Table. Values of the functions ej(to), fiito), cii(to), hito), Cj(to), and di(to) versus a 

choice of x^(0), A^(0), and ^^(0). 


x^(0) 

0 

0 

0 

0 

0 

4 

10 

10 

A'^(O) 

2 

4 

10 

10 

10 

4 

10 

10 

h'm 

10 

10 

10 

4 

2 

10 

10 

4 

et{to) 

0.0011 

0.0024 

0.0082 

0.0310 

0.0651 

0.0070 

0.0102 

0.0308 

ex{to) 

0.2127 

0.1192 

0.0260 

-0.0227 

-0.0544 

0.0661 

0.0113 

-0.0232 


0.5066 

0.3451 

0.1804 

0.0774 

-0.0062 

0.0430 

-0.0136 

-0.0528 

ft{io) 

1.7196 

1.7406 

1.7698 

1.8352 

1.9272 

1.7254 

1.7489 

1.8119 

/a( to) 

-0.4414 

-0.4295 

-0.4161 

-0.3976 

-0.3722 

-0.4386 

-0.4262 

-0.4046 

f>c{to) 

0.0173 

0.0195 

0.0175 

-0.0069 

-0.0422 

0.0232 

0.0273 

0.0138 

ai(to) 

0.7533 

0.6725 

0.5902 

0.5387 

0.4969 

0.7632 

0.6882 

0.6247 

0-2{to) 

-0.6217 

-0.6601 

-0.6927 

-0.6842 

-0.6539 

-0.6079 

-0.6406 

-0.6415 

(ly{to) 

0.5066 

0.3451 

0.1804 

0.0774 

-0.0062 

0.0430 

-0.0136 

-0.0528 

OQ(to) 

0.5417 

0.5558 

0.5724 

0.5924 

0.6164 

0.5430 

0.5571 

0.5779 

auito) 

0.0833 

0.1116 

0.1448 

0.1847 

0.2328 

0.0859 

0.1141 

0.1558 

bi{to) 

0.5557 

0.5665 

0.5761 

0.5757 

0.5703 

0.5538 

0.5627 

0.5646 

h(to) 

-3.0399 

-3.0070 

-2.9664 

-2.9044 

-2.8175 

-3.0352 

-3.0014 

-2.9402 

hyito) 

0.0724 

0.0939 

0.1131 

0.1122 

0.1015 

0.0695 

0.0907 

0.1004 

bgito) 

5.3129 

5.3202 

5.3305 

5.3514 

5.3821 

5.3150 

5.3234 

5.3432 

bu{to) 

3.6951 

3.7099 

3.7305 

3.7722 

3.8336 

3.6995 

3.7162 

3.7557 

Cl(to) 

0.0034 

0.0036 

0.0040 

0.0087 

0.0145 

0.0024 

0.0034 

0.0080 

C2(to) 

-0.0085 

-0.0096 

-0.0165 

-0.0543 

-0.1080 

-0.0159 

-0.0197 

-0.0549 

Cj/(to) 

0.0068 

0.0072 

0.0079 

0.0173 

0.0290 

0.0052 

0.0056 

0.0127 

CQ(to) 

-0.0040 

-0.0044 

-0.0068 

-0.0210 

-0.0409 

-0.0061 

-0.0077 

-0.0210 

Cf/(to) 

-0.0079 

-0.0088 

-0.0137 

-0.0420 

-0.0817 

-0.0122 

-0.0154 

-0.0419 

di(to) 

-0.0186 

-0.0144 

-0.0055 

-0.0008 

0.0023 

-0.0063 

-0.0029 

-0.0010 

d2(to) 

-0.0143 

-0.0114 

-0.0069 

-0.0122 

-0.0217 

-0.0076 

-0.0062 

-0.0129 

dyi^to^ 

-0.0372 

-0.0288 

-0.0109 

-0.0016 

0.0045 

-0.0196 

-0.0034 

0.0044 

dqito) 

0.0014 

0.0010 

-0.0005 

-0.0038 

-0.0080 

-0.0004 

-0.0011 

-0.0039 

duito) 

0.0029 

0.0020 

-0.0010 

-0.0076 

-0.0160 

-0.0008 

-0.0022 

-0.0079 

(^o) 

0.0033 

0.0073 

0.0245 

0.0929 

0.1953 

0.0210 

0.0305 

0.0923 

hA(to) 

0.6382 

0.3575 

0.0779 

-0.0680 

-0.1631 

0.1983 

0.0340 

-0.0695 


1.5199 

1.0352 

0.5411 

0.2323 

-0.0185 

0.1290 

-0.0407 

-0.1583 

bt{to) 

5.9680 

6.0231 

6.0947 

6.2192 

6.3981 

5.9794 

6.0382 

6.1588 

&A(to) 

-2.4118 

-2.3466 

-2.2771 

-2.2165 

-2.1457 

-2.4119 

-2.3479 

-2.2752 

bU^o) 

0.2172 

0.2817 

0.3394 

0.3367 

0.3045 

0.2085 

0.2722 

0.3011 

Cl (to) 

-0.0204 

-0.0228 

-0.0370 

-0.1173 

-0.2306 

-0.0342 

-0.0427 

-0.1178 

CA(to) 

0.0017 

0.0012 

-0.0047 

-0.0283 

-0.0645 

-0.0083 

-0.0107 

-0.0342 

(^o) 

0.0204 

0.0216 

0.0238 

0.0520 

0.0870 

0.0156 

0.0166 

0.0380 

(^o) 

-0.0099 

-0.0084 

-0.0084 

-0.0235 

0.0457 

-0.0088 

-0.0095 

-0.0247 

dxito) 

-0.0700 

-0.0547 

-0.0233 

-0.0146 

-0.0150 

-0.0334 

-0.0125 

-0.0096 

dy^(to) 

-0.1115 

-0.0865 

-0.0327 

-0.0049 

0.0135 

-0.0587 

-0.0103 

0.0132 
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Figure captions 


Fig. 1. Evolution of the trilinear couplings At{t) and Ax{t) from the Grand Unihcation 
scale (t = 0) to the electroweak scale (t = to) at = 0 and A^(0) = h^{0) = 10 for 
universal initial conditions. Parameter A varies in the range —Mi 12 < A < M1/2. 


Fig. 2. Evolution of the combinations of masses 911? (t) and 911? (t) from the Grand Uni¬ 
hcation scale {t = 0) to the electroweak scale {t = to) at = 0 and A^(0) = h?(0) = 10 
for universal initial conditions. Parameter m? varies in the range 0 < m? < 

Fig. 3. Evolution of the trilinear couplings At(t), Ax(t), and from the Grand 

Unihcation scale (t = 0) to the electroweak scale (t = to) at h?(0) = A^(0) = x^(0) = 10 
for universal initial conditions. Parameter A varies in the range —M 1/2 ^ A < Mi/ 2 - 

Fig. 4. Evolution of the combinations of masses 911? (f), 911?(f), and 91t?.(t) 
from the Grand Unihcation scale (t = 0) to the electroweak scale (t = to) at 
ft|(0) = \H0) = x^(0) = 10 for universal initial conditions. Parameter m? varies in the 
range 0 < m? < 

Fig. 5. Boundary conditions for the renormalization group equations of the NMSSM 
at the Grand Unihcation scale {t = 0) at x^ = 0 and h?(0) = A^(0) = 20 uniformly 
distributed in the {At/ Mi 121 Ax/ Mi 12 ) plane - Fig. 5a, and the corresponding values of 
the trilinear couplings at the electroweak scale {t = to) - Fig. 5b. The straight line in 
Fig. 5b is a ht of the values (At(fo), ^A(to))- 


Fig. 6 . Boundary conditions for the renormalization group equations of the NMSSM at 
the Grand Unihcation scale (t = 0) at x^ = 0, hf(0) = A^(0) = 20, and At(0) = Aa(0) = 0 
uniformly distributed in the (DJtf/M/y, DJlf/M/y) plane - Fig. 6a, and the corresponding 
values of the trilinear couplings at the electroweak scale (t = to) - Fig. 6b. The straight 
line in Fig. 6b is a ht of the values (9Jl?(to), 91l?(to)). 

Fig. 7 . Planes in the parameter spaces (At/Mi/2, Aa/Mi/ 2, A^/Mi/2) - Fig. 7a., and 
(W^t/M/y, W^l/M/y, V!}ll./M/y) - Fig. 7b. The shaded parts of the planes correspond 
to the regions near which the solutions at h?(0) = 16, A^(0) = 12, and x^(0) = 6 are 
concentrated. The initial values Aj(0) and 91l?(0) vary in the ranges —M 1/2 < A < M 1/2 
and 0 < 911?(0) < 3M/y, respectively. 

Fig. 8 . Planes in the parameter spaces (At/Mi/2, Ax/Mi/2, A^/Mi/2) - Fig. 8a, and 
(fXnf/M/y, W’ll/M/y, V7ll./M/y) - Fig. 8b. The shaded parts of the planes correspond 
to the regions near which the solutions at h?(0) = 32, A^(0) = 24, and x^(0) = 12 are 
concentrated. The initial values Aj(0) and 91l?(0) vary in the ranges —M 1/2 < A < M 1/2 
and 0 < 91l?(0) < 3M/y, respectively. 
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Fig. 9. Set of points in planes (0.0223(A^/Mi/2)+0.1278 (Aa/Mi/2), A/^ 1 / 2 ) ^ Fig. 9a, 
and (0.0223(97l^/M^y2) + 0-1278(07l^/M^^^2)) ^t/^ 1 / 2 ) “ corresponding to the 

values of parameters of soft SUSY breaking for h‘f{0) = 32, A^(0) = 24, x^(0) = 12, 
and for a uniform distribution of the boundary conditions in the parameter spaces 
{At, Ax, and (971^, The initial values A(0) and 971^(0) vary in the ranges 

—Mi /2 < ^ < Mi /2 and 0 < 971^(0) < ‘iM’li^, respectively. The straight lines in Figs. 9a 
and 9b correspond to the planes in Figs. 8a and 8b, respectively. 


26 



^x/Mi/ 2 ^ At/M i/2 



Fig.la. 



Fig.lb. 
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